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Abstract. The CGMV method allows for the general discussion of localization 
properties for the states of a one-dimensional quantum walk, both in the case 
of the integers and in the case of the non negative integers. Using this method 
we classify, according to such localization properties, all the quantum walks with 
one defect at the origin, providing explicit expressions for the asymptotic return 
probabilities to the origin. 



1. Introduction 

A quantum random walk can be considered as a quantum analog of the more 
familiar classical random walk on a lattice. In this much simpler case the study of 
interesting "return properties" can be said to have started with G. Plya (1921), |31j . 
who proved that the simplest unbiased walk will eventually return to the origin with 
probability one in dimension not greater than two. This holds in spite of the fact 
that the probability of returning to the origin in n steps, denoted by p(n), converges 
to zero as n tends to infinity. 

In this paper we consider aspects of this problem in the context of quantum 
random walks (QWs). We give a method that allows us to analyze the asymptotic 
behaviour of the quantity p(n) for two-state one-dimensional QWs, leading to the 
discovery of general features of this asymptotics in the case of distinct coins. We 
also apply this method to the QWs that are given by one arbitrary common coin 
at each site except for an arbitrary "defective" coin at the origin. Before giving a 
summary of the results in the paper we give a brief review of the more standard 
case of a classical random walk. 

In this more traditional case, and for the simplest unbiased walk, one can obtain 
an expression for p(n) in many different ways. This is basically true since one is 
dealing with a translation invariant evolution. As soon as this condition is relaxed 
things become much harder. One of the methods that can (at least in theory) give an 
expression for p(n) goes back to the work of S. Karlin and J. McGregor (1959), [23] . 
Their method applies to a birth-and-death process on the nonnegative integers but 
they themselves already contemplated extending their method to such processes on 
the integers by using matrix valued objects. This has been implemented recently 
and independently by H. Dette et al. (2006), [9], as well as by FA. Griinbaum 
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(2007), [12]. In both cases one makes crucial use of the matrix valued orthogonal 
polynomials introduced by M.G. Krein (1949), [281125]. 

The Karlin-McGregor (KMcG) method alluded to above proceeds by starting 
from a very simple case, say a situation with no defects, considering its orthogonality 
measure and related function theoretical objects such as its Stieltjes transform. One 
then introduces one defect and studies the effect that this has on the orthogonality 
measure. In this way one can get an expression for the new value of p(n). Clearly 
this process can be iterated a finite number of times to obtain situations that are 
far from the initial basic case. The KMcG method proceeds by studying the effect 
of the defect on the Stieltjes transform of the measure and it produces much more 
than an expression for the new value of p(n). The interested reader can find a host 
of examples treated in the fashion described above in the papers [HI M 021 E2 EH 

uni [it] . 

In a previous paper, [5], we found an analog of the KMcG method that can 
be used in the case of QWs on the integers and the non negative integers. This 
method has recently been used by other authors, see [27] . We use their terminology 
and call this the CGMV method. In the CGMV method the central role is played 
by a convenient spectral representation of the unitary operator that describes the 
dynamics of the walk. In that sense this is very close in spirit to the mathematical 
foundations of Quantum Mechanics laid down by people like E. Schredinger, W. 
Heisenberg, M. Born and J. von Neumann. 

The connection with QWs is embodied in the fact that a unitary operator has a 
nice five-diagonal representation, see [I]. While this is true for any unitary operator, 
the search for such a five-diagonal representation can be a difficult task in the infinite 
dimensional case. However, for a standard two-state one-dimensional QW the five- 
diagonal representation comes from an appropriate reordering of the usual basis of 
states. 

While the method itself is fairly similar to the one used by Karlin and McGregor, 
the results that one obtains in the quantum case are, regardless of the method, often 
intrinsically different from the classical ones. For instance, the application of CGMV 
tools yields in [5] a number of situations dealing with the notion of recurrence, where 
one sees that classical and quantum walks have little in common. Some of these 
differences were previously discovered in [351 ESI |37] using a Fourier approach to 
quantum recurrence. 

It is possible that many more such differences will become apparent as one devel- 
ops good tools to study these kinds of walks. For instance, appropriate mathematical 
techniques such as the fractional moment method, [TJ, have been key to prove for 
QWs with random coins a peculiarity of quantum systems in random environments 
called "Anderson localization": the wave packets stay trapped in a finite region of 
space for all time. There is a huge literature on this subject that has its origin in a 
discrete model in solid state physics, see [3] . The interested reader may get a guide 
to the literature as well as very nice discussion of this notion by consulting [18] . 
Concerning Anderson localization in QWs see [22] . 

For an environment that is not too disordered, weaker "localization" properties 
could distinguish quantum from classical random walks too. A candidate for such 
a property is related to the asymptotics of the return probability to a given site. 
To be more precise we start at a qubit state (a, (3) at a site k on the lattice and 
look at the probability p^\(n) that after time n the state will again be localized 
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at site k but with unspecified spin orientation. We adopt the terminology of [26] 
and say that the qubit (a, 0) at a site k exhibits localization when p^\{n) does not 
converge to zero as n tends to infinity. We will deal with this notion of localization 
of single qubit states, which should not be confused with such a global property 
of disordered quantum systems like Anderson localization, and sometimes we will 
refer to it as "single state localization" . 

Irreducible two-state QWs on the integers with a constant coin can not exhibit 
localization, see [2] for instance. Nevertheless, localization can appear for homoge- 
neous QWs if we increase the internal degrees of freedom, [201 HH] , or the dimension 
of the lattice, [2~T| l38j. A way to get localization in one dimension keeping a constant 
coin is to destroy the translation invariance by considering the lattice of the non 
negative integers, see [51 127]. 

Another way of breaking the translation invariance is to introduce defects on the 
integers. In this context an important difference between classical and quantum 
random walks has already been recognized concerning localization of single states. 
In [26], N. Konno has studied a perturbation of the Hadamard QW on the integers 
and proved that in the case of a particular defect at the origin, and starting from 
a particular qubit too, the quantity p^\{ n ) ^ a ^ s t° converge to zero. The method 
used in [26J is a path counting argument. One can think of our paper as an effort 
to explore the phenomenon uncovered in [26] in a more general setting by using the 
CGMV method. 

Localization can depend on the initial state as well as on the coins of the QW. 
Our aim is to get a better understanding of these dependencies. This appears to be 
a difficult task if one resorts to the usual approaches, like path counting or Fourier 
transform, which usually allow for the computation of the asymptotics of p^ain) at 
specific initial qubits or for very limited coin models. 

The state and coin dependencies of localization seem to be more tractable from 
the CGMV point of view. This is specially true for the state dependence, and we 
will be able to have a picture of it for a large class of QWs. The coin dependence, 
which is much more involved, will be discussed in some explicit examples. More 
precisely, the CGMV approach will allow us to perform a complete classification, 
according to the localization behaviour, of all the QWs with a coin which is constant 
except for a defect at the origin. 

In addressing the study of the dynamics of any quantum system one should keep 
in mind that there is a large literature on the subject. The classical book "Non- 
relativistic quantum dynamics" by W. Amrein (1981) gives a good treatment. For 
a more recent account see the book "Hilbert space methods in quantum mechanics" 
(2009) by the same author. 

Most of the work deals with the spectral properties of the unitary group that 
implements the quantum evolution and its dynamical consequences. There are also 
several papers dealing with these issues of which we just mention two: at a fairly 
technical level one can consult [30] . or later work of this same author, and at a more 
approachable level see [22]. 

In a number of ways one can say that what we do in great detail is related to 
the bread and butter of Quantum Mechanics. What we call the return probabil- 
ity p^\{ n ) is related (but not identical) to what Y. Last, [30], calls the "survival 
probability", see (1.3) of his paper. The main point of our paper is that we manage 
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to compute these quantities explicitly and then study their asymptotic values as n 
goes to infinity 

We now try to give an account of the contents of the present paper, which deals 
with the localization properties of two-state QWs on the integers and the non neg- 
ative integers. 

Just as in the KMcG method used for classical random walks one studies the 
effect that introducing defects on a simpler walk has on the Stieltjes transform of 
the orthogonality measure, in the quantum case we need to study the so called 
Schur function of the measure. The analysis of the special features of the Schur 
functions related to QWs with distinct coins, with special emphasis in the case of 
the integers, takes up a good part of section [2] in the present paper. The results of 
this section will be key for the rest of the paper. 

The general problem of the single state localization within the CGMV language is 
studied in section |3j The main result is Theorem [331 which establishes a connection 
between localization in a QW and the singular part of the corresponding orthog- 
onality measure. In particular, the absence of such a singular part leads to the 
absence of localized states, while the presence of mass points implies the existence 
of states which exhibit localization. 

Among other consequences the CGMV method shows that, despite its name, sin- 
gle state localization is in fact a quasi global property for a large class of QWs so 
that a localization dichotomy holds in many cases: either no state exhibits localiza- 
tion or at most one state per site is localization free. Such a dichotomy is ensured 
when the singular part of the measure is not purely continuous. That is, the state 
dependence of localization is quite regular for a wide class of QWs. In this case we 
will refer to QWs with or without localization omitting any mention to the initial 
qubit state. 

In section [3] we see that the localization dichotomy holds in particular for any QW 
with periodic coins up to a finite number of defective coins. It is also shown that, 
among these QWs, the case of strictly periodic coins on the integers is somewhat 
special because it never gives localization. QWs on the integers whose coins have 
period P are the only one-dimensional QWs which are invariant with respect to 
right and left translations of P sites. Thus, we can state that the requirement of a 
(right and left) translation invariance for a one-dimensional QW forces the absence 
of localization. 

In sections H] to [8] these results are made more specific, both for the non negative 
integers as well as for the integers, in the case of the simplest perturbation of the 
constant coin model: the QWs with a coin which is constant except for one defect 
at the origin. These QWs will serve as a laboratory to study the coin dependence 
of localization in QWs. 

As we pointed out, localization already appears for two-state QWs with a constant 
coin on the non negative integers but not on the integers where the related measure 
is absolutely continuous (see [5]). Therefore, the study of QWs with one defect 
acquires a special relevance in the case of the integers because they are the nicest 
laboratory in which one can study localization on such a lattice. Nevertheless, we 
will perform also the analysis of one defect on the non negative integers, which will 
allow us to compare with the case of the integers, thus showing the effects of the 
boundary conditions on the localization behaviour. 
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As a particular case of the periodic QWs with a finite number of defects, the 
localization dichotomy holds for QWs with one defect. The analysis of localization in 
these models becomes the study of the presence of mass points in the corresponding 
measure. 

Section [4] deals with the general features of the orthogonality measure for QWs 
with one defect at the origin. It is shown that they fall into groups with the same 
measure up to rotations, thus with the same localization behaviour. These groups 
are labelled by two parameters a, b in the unit disk for the non negative integers, 
while an additional labelling parameter u in the unit circle appears in the case of 
the integers. 

Section shows that u actually plays no role in the presence or absence of localiza- 
tion, hence localization for one defect at the origin only depends on two parameters 
a, b which are defined by the coins of the QW (in a different fashion for the integers 
and the non negative integers, see (THj) and ([191) ). The parameter a depends only 
on the unperturbed coin and the phases of the perturbed one, while b depends on 
the perturbed coin and the phases of the unperturbed one. 

Sections [5] and [6] give a very exhaustive analysis of localization for one defect on 
the integers, and the same in depth analysis is carried out in sections [7J and [8] for 
the case of the non negative integers. Sections \5\ and [7J discuss the coin dependence, 
providing a characterization of localization in terms of the parameters a, b. Sections 
M and [S] yield explicit results for the asymptotic return probability to the origin for 
any defect and any initial qubit state. 

Different localization figures in the space of parameters a, b are presented in sec- 
tions [5] and [7J They demonstrate that, in contrast to the classical case (see [261 
Section 6]), localization is dominant under the presence of a defect. Nevertheless, 
these figures also show that, at the same time, there are situations where the ab- 
sence of localization is stable under small perturbations of a and b, i.e, under small 
perturbations of the coins. In particular, given \a\, the largest regions for the pa- 
rameter b without localization appear when a is imaginary, both for the integers 
and the non negative integers. Then there is no localization if Imfe > Ima > or 
Imfe < Ima < 0. 

Concerning the return probabilities p^\(n), the CGMV method not only shows 
the dependence of its asymptotics on the initial qubit (a,j3), but also explains the 
reason for its possible oscillatory asymptotic behaviour: the presence of the factors 
z n in (fTTj) . where z are the mass points of the measure. The return probabilities 
turn out to be convergent when the singular part of the measure is a unique mass 
point or, in the case of several mass points, when some symmetries force the mutual 
cancellation of the cross terms in (Hip . 

General reasons imply that the return probabilities p^L(2n— 1) at odd time vanish 
for any QW on the integers, which is related to the fact that the mass points on 
the integers always appear in pairs which are symmetric with respect to the origin. 
Therefore, in the presence of localization on the integers we can not expect the 
convergence of Pas( n ) but at most of p^(2n). This convergence takes place for 
sure when the singular part of the measure is a single pair of symmetric mass points. 

QWs on the integers with one defect at the origin have a symmetry under reflec- 
tion of the sites with respect to the origin which causes the convergence of p^\(2n) 



6 



CGMV 



regardless of the number of mass points. However, for one defect on the non nega- 
tive integers with more than one mass point, as well as when considering the return 
probability p^(2n) to a site k ^ on the integers with more than two mass points, 
the asymptotic behaviour is in general oscillatory. 

The state dependence can also disappear in some special situations like, for in- 
stance, one defect at the origin on the integers with an imaginary value of a. In 
this case section [6] proves that p^p(2n) actually converges to the same limit for any 
initial qubit. This covers as a special case the result obtained in [26] for a concrete 
perturbation of the Hadamard model and a specific initial state. We not only prove 
that the result in [26] is state independent, but we also extend this to a more general 
model with one defect since we find that the state independence holds whenever the 
products of the diagonal coefficients of the perturbed and the unperturbed coins 
have the same phase. 

Section [6] gives explicitly p^ = lim n _ s . 00 p^(2n) for any QW on the integers 

with one defect at the origin. The convergence of p^(2n) allows for the analysis 
of the maximum asymptotic return probabilities to the origin, both when running 
over the qubits (a,/3) and also when running over the parameters a,b, i.e., over 
the coins of the model. We find that max a) pp^\ approaches one when \a\ — > 1 
provided that |Ima — Im6| is bounded from below. The consequence is that, given 
a defective coin, for most of the choices of the non defective coin there exist qubits 
which asymptotically return to the origin with probability almost one, as long as 
the non defective coin is close enough to an anti-diagonal one. 

In the special case of an imaginary a, the asymptotic return probability p ®n does 
not depend on the state and approaches one when \a\ — > 1 if \a — b\ is bounded from 
below. This implies that, when the products of the diagonal coefficients of both 
coins have similar phases, all the qubits asymptotically return to the origin with 
probability almost one, provided that the non defective coin is close enough to an 
anti-diagonal one. 

The results above not only show the strength of the CGMV method for the anal- 
ysis of localization in QWs, but they address new research lines with could lead to 
new and surprising quantum effects. For instance, it could be very interesting to 
analyze the localization behaviour of QWs where the localization dichotomy is not 
ensured, i.e., those whose measure has a singular part which is strictly continuous, 
and specially those with a purely singular continuous measure. The physical conse- 
quences of singular continuous spectra in Quantum Mechanics is an active field of 
research, see for instance [30] and the references therein. The study of this problem 
in those models which can be considered the simplest realization of a dynamical 
quantum system, i.e., the QWs on a lattice, could make it easier to understand the 
quantum meaning of a singular continuous spectrum and its dynamical implications. 



Throughout the paper we will deal with QWs on a state space with a basis 
{|fc t), \k i)}kez or z+j where Z + = {0, 1, 2, . . . }. The quantum dynamics will be 
governed by unitary coins 



2. QWS, CMV MATRICES AND SCHUR FUNCTIONS 
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so that the one-step transitions are given by the operator il defined by 

U|fcT) = cS fc 1 ) |A; + lt) + 4 fe 1 ) |A;-H), U|*4) = c$\k + It) + c$\k - U), 
except when = for Z + , in which case the unitarity forces a transition 

it|0t> = cg>|lt> + 4°i ) |0t), U|(4) = cSllt) + 4° 2 ) |0t). 

If a diagonal element of a coin were null, the QW would decouple into independent 
ones, so the requirement ^ is not really a restriction, but it simply means 
that we are considering only irreducible QWs. 

Once an order is chosen for the basis {|fc t), \ k 4)}fceZorZ + , it gives a matrix 
representation U of the transition operator il, which we will call the transition 
matrix of the QW. If \z) denotes the i-th vector of such an ordered basis, we will use 
the convention that U = {U^j) is defined by = J2j so that the one-step 

evolution — > reads as if) — > if)U using the coordinates if) = (if)o, if>i, . . .) of 

l*) = £^)- 

It was shown in [5] that the order 
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Afc, CKfc and p£'" being scalars for Z + and 2 x 2-matrices for Z. 

If we refer at the same time to both, scalar and matrix objects, we will use the 
boldface notation. However, whenever we wish to distinguish between scalars and 
matrices we will reserve for the last ones the boldface notation, so in such a case we 
will denote by \ k , a k and Pk = P k = Pk the scalars of a QW on Z + . 

. (fc) (k*\ 

Denoting by e lf7 j the phase of , the coefficients X k , cx k and p k are obtained 
from the coins by means of 
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The matrix C is a special case of a kind of matrices that play a central role in the 
theory of orthogonal polynomials (OP) on the unit circle T = {z G C : \z\ = 1} , the 
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so called CMV matrices 
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The case of a QW corresponds to 0:2/0+1 = so that p^+i — P2I+1 = 1 stands for 
the number 1 or the 2x2 identity matrix for Z + and Z respectively. 

The connection between the transition matrix U of a QW and the CMV matrices 
ensures that the Laurent polynomials X k defined by 

UX(z) = zX(z), X = (X ,X 1 ,...) T , X (z) = l, (3) 

constitute a sequence of orthogonal Laurent polynomials (OLP) with respect to a 
probability measure /x supported on the unit circle (see [5]), i.e., 

/ X,i(z)dii(z)X k (zy = 15 jjk . 
Jt 

Such OLP and measures are scalar or 2 x 2-matrix valued for Z + and Z respectively. 
The coefficients a. k are known as the Verblunsky or reflection coefficients of the 
measure [i. 

The orthogonality and the "eigenvalue" equation ([3]) yield a KMcG formula for 
the QW, i.e., an OLP representation of the ra-step transition amplitudes (see [5], 
pages 479 and 483]) 

(U n ) hk = [ z n Xi{z)dn{z)X k {z)l 
Jt 

Here (•) Ji fc stands for the (j, k)-th element in Z + and for the (j, k)-th 2 x 2-block 
in 7L. The KMcG formula is the cornerstone of the CGMV method, which takes 
advantage of the OP techniques for the analysis of QWs. 

Useful tools in the theory of OP on the unit circle are the so called Carathodory 
and Schur functions related to ^t, defined respectively by 

F(z)= [t±±drtt), f{z) = z- 1 {F{z)-l){F{z) + l)'\ \z\<\. 

The Carathodory and Schur functions of a probability measure on T can be charac- 
terized as the analytic functions on the unit disk = {zGC:|,2:|<l} such 
that F(0) = 1, ReF(z) > and \\f(z)\\ < 1 for z G D respectively, where 
KeA — 2 (A + A^) and ImA = ^(A — A*) for any square matrix A. We will assume 
that F and / are extended to the unit circle by their radial boundary values, which 
exist for almost every point of T, so ReF(e l6 ) > and ||/(e j6l )|| < 1 for a.e. 9. 

These analytic functions provide a shortcut between the measure /j, and the 
Verblunsky coefficients cx k because both of them can be recovered from F or /. 
We will now point out the connections needed for the rest of the paper (see [S] for 
the general matrix case). 



ONE-DIMENSIONAL QUANTUM WALKS WITH ONE DEFECT 9 

If dp(e %e ) = w(9)t^ + dp s (e id ) is the Lebesgue decomposition of p into an abso- 
lutely continuous and a singular part, 

w{9) = ReF(e w ), for a.e. 9, 
det w(6) ^ \\f(e ie )\\ < 1, for a.e. 9, 
supp/x s C {z G T : linvri tr(ReF(rz)) = oo}, 
p({z}) = lim rtl ^F(rz), zGl 

In particular, if z G T is a pole of an analytic extension of F, such a pole must be of 
order one and z is an isolated mass point of p with mass p({z}) = — (2z)~ l ~Res(F; z). 

These properties are well known in the case of scalar measures. Concerning 
matrix measures, the first two properties can be found in [8], while the remaining 
ones can be reduced to the scalar case by noticing that p is absolutely continuous 
with respect to the scalar trace measure tip. 

On the other hand, starting at / = /, the Verblunsky coefficients oc k = f k (0) 
can be recovered through the Schur algorithm 

f k+1 (z) = z-\p*)-\f k {z) - <x k ){l - cJj k {z)y l p L k 

= z-'p* (1 - f k {z)oc\)-\f k {z) - a h ) (p L k r\ 

which assigns to / a sequence of Schur functions f k . For this reason, the Verblunsky 
coefficients of p are also known as the Schur parameters of /. Obviously, the Schur 
parameters of each Schur iterate f k are obtained deleting the first k parameters of 
/: ct k , ct k+ i, .... The Schur algorithm can be inverted to give 

f k ( Z ) = (p«r\zf k+1 (z) + a fc )(i + a \zf k+1 {z))- i P L k 

= pi (1 + zf k+1 (z)<x{)-\zf k+1 (z) + a*) (ft)- 1 . 

In the scalar case the factors p k cancel each other in the above formulas due to 
the commutativity, giving 

1 f k {z)-a k zf k+1 (z)+a k 

f k+ i[z) = — _ , fk{z) = t-t. 5 

zl-a k f k (z) l + a k zf k+ i(z) 

Both, the measure and the Schur function, are univocally determined by the 
Verblunsky coefficients. Some results relating Schur functions and Schur parameters 
will be of interest for us. 

Proposition 2.1. If cto, a.2, cx^, . . . are the Schur parameters of the Schur function 
f(z), then et , 0, ck 2 , 0, a 4 , 0, . . . are the Schur parameters of f(z 2 ). 

Proof. Denote g 2k -i( z ) = z fk( z2 ) an d 92k( z ) = fk( z2 )- The substitution z — > z 2 
in the Schur algorithm for f(z), together with the relation g 2k ( z ) = z ~ 1 92k-i( z )i 
gives the Schur algorithm for g(z) = f(z 2 ). □ 

Proposition 2.2. If a Schur function f has Schur parameters cx k , then for any 
unitary matrices Vi, V 2 the Schur function V\fV 2 has Schur parameters ViCt k V2- 

Proof. Simply check that the Schur algorithm is invariant under the transformation 
/ — > VifV 2 , cx k — > ViCt k V 2 , bearing in mind that it maps p k — > V 2 p k V 2 and 

p* -> v lP *vl □ 
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Proposition 12.11 states that the Schur functions whose odd Schur parameters van- 
ish are the even Schur functions, i.e., those satisfying f(—z) = f(z). In terms of 
the Carathodory function this condition reads as F(—z)F(z) = 1, as follows from 
the inverse relation 

F(z) = (l + zf(z))(l-zf(z))- 1 . (6) 

Hence, the Schur function of any QW on Z or Z + must be even. On the other 
hand, Proposition 12.21 has the following consequences of interest for QWs on Z. 



a k 
at 



Proposition 2.3. Given a sequence of 2x2 Schur parameters 
the corresponding Schur and Carathodory functions are 

™ " (u% '-o W 

where f± is the Schur function with Schur parameters ■ 

If dfi(e ie ) = w(9)y- + d[i s (e ie ) is the Lebesgue decomposition of the related mea- 
sure, 

detw(e id ) ^ \f±(e ie )\ < 1, for a.e. 9, 

and the singular part n s is supported on the roots z G T of g(z) = 1. The mass 
points are those roots such that 

1 — r 

m(z) = lim ; — r ^ 0, 

rti 1 - g(rz) 

and the corresponding mass is the singular matrix 

/*({*}) = m(z) (-Lj , V (z) = zf.(z) G T. 

In particular, if g extends analytically to a neighbourhood of a root z G T of g(z) = 
I, then z is a simple isolated roo£\ and also an isolated mass point with m(z) = 
l/zg'{z). 

Proof. The result for / is a direct consequence of Proposition 12.21 and 

f+ o\ _ o \ „_ fo 1 



f = v^ f _y o, = v^ ' a -y v- v 

Then, the expression of follows from ([6]). 

The rest of the results are obtained from (j3J) by using the actual form of the Schur 
and Carathodory functions. Concerning the singularity of the masses, simply take 
into account that the roots z G T of g(z) = 1 must satisfy |/±(z)| = 1 because 
|/±| < 1 in T. Hence, the non diagonal elements of the mass must be proportional 
to 7/0) = zf-(z) G T and zf+(z) = (zf^(z))- 1 = : q{z). □ 



1 If the root z were not isolated, general principles would imply that g = 1 on D, which is 
not possible since g must be a Schur function. Similar comments for the roots of h = 1 in the 
paragraph below Corollary 12.41 
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Proposition 12.31 holds for QWs on Z with a\ k = a 2 k, a 2k = —a_ 2 k-2 and &tk+i = 

a 2k+i = 0' where a 2 k = c^i ^k/Mk-i- Then, / + is the Schur function associated 
with the Schur parameters with non negative indices, while /_ corresponds to the 
Schur parameters with negative indices. Furthermore, f± are even functions, so g 
is even too, and this has the following consequence. 

Corollary 2.4. For any QW on Z, the mass points of the corresponding measure 
appear in pairs ±z which are symmetric with respect to the origin, and the mass is 
given by Proposition \2.3\ with m(—z) = m(z) and T)(—z) = —r)(z). 

The mass points of a QW on Z + have no such a symmetry, despite the fact that 
the corresponding Schur function / is even too. The reason is that the Carathodory 
function is given by 

F{z) = h{z) = zf(z), 

so the singular part of the measure is supported on the roots z G T of the equation 
h(z) = 1, which is not invariant under z — > —z because h is odd. The mass points 
are those roots such that 

1 — r 

li({z}) = lim — — r ^ 0. 

rti 1 - h{rz) 

When h has an analytic extension to a neighbourhood of a root z £ T of h(z) — 1, 
such a root is simple and isolated, and is also an isolated mass point with mass 
»({z}) = l/zh'(z). 

3. Single state localization in QWs 

Following [26], we will adopt the definition below for the localization of a state 
in a QW. It applies only to qubit states a\k t) + 0\k \) at a given site k, and 
characterizes those states which have a non null probability of asymptotic return to 
the same site where they are placed originally. 

Definition 3.1. Given a QW on Z or Z + , let p^g(n) be the probability that the 

walker returns to the site k in n steps having started at the qubit state l^^jg) = 
a\k^\) + f3\k \) at the initial time. We will say that such a state exhibits localization 
if limsupn^oopj^n) ^ 0. 

It is known that the structure of the transition matrix for the QWs on Z we 
are discussing always gives a null return probability for an odd number of steps, 
i.e., p^\{2n — 1) = 0. Therefore, the only quantity of interest in the case Z is the 

asymptotics of p^(2n). 

If ii is the transition operator of the QW, 

The KMcG formula provides an alternative expression for this probability which is 
nicely adapted to study its asymptotics. Indeed, a simple extension of the formula 
in page 497] to the case of two arbitrary states gives 

(*|lt n |^) = i)TJ n ^ = [ z n ^(z)d f i(z)^(z)\ (7) 

Jt 
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where if>(z) is an L 2 (Y) function associated with the state = Ylii^i^) (K) * s the 
i-th vector of the ordered basis), which is a scalar function for Z + and a 2- vector 
function for Z. The general form of if)(z) is given in the first column of the following 



table, while the second column shows the particular case t/>?g(z) for the qubit state 

i*S>. 









Z+ 


k 


aX 2k + (3X 2 k+i 


z 


^2(ljj2k,1p2k+l)X k 

k 


\ 


\(a,0)X 2j + (0,p)X 2j+1 k = j 
(0, (3)X 2j + (a, 0)X 2j+1 k = -j - 1 3 ~ 



Relation d7J gives the identity 



to 




+ 


/ 



z^{z)d^(z)^>(z) 



The importance of the KMcG formula in the study of the localization of the states 
in a QW was first pointed out by N. Konno et al in [27]. There the authors use 
the Riemann-Lebesgue lemma to obtain the asymptotics of p^{n) for the case of 
a constant coin on Z + . The method can handle other QWs on Z + , as well as QWs 
on Z. 

For convenience, in what follows we will use the notation 

a n ~ K lim (a n - b n ) = 0. 

n n— >oo 

Lemma 3.2. If 51 is the transition operator of a QW on Z or Z + with measure [i, 



n J j 



where fi s is the singular part of fi. 

Proof Let dfj,(e ie ) = w(6)^ + dfj, s (e ie ) be the Lebesgue decomposition of the mea- 
sure \x. The Riemann-Lebesgue lemma implies that 



lim 



27T _ in 

e in6 ^{e ld )w{e)^{e ie ) ] — = 

2ir 



because %jj{e )w(9)tp(e y is integrable with respect to the Lebesgue measure. This 
gives the result, bearing in mind the KMcG formula. □ 

As a consequence, no state can exhibit localization in a QW with an absolutely 
continuous measure. As for the singular part, it always can be decomposed into mass 
points and a singular continuous part. As we will see, due to Wiener's theorem, the 
presence of mass points will always give localized states, regardless of the presence 
of a singular continuous part. However, if the singular part is exclusively continuous 
the situation is more involved because the Riemann-Lebesgue lemma holds for some 
singular continuous measures, but not for all of them. 

To obtain the strongest results about localization for QWs on Z one is greatly 
aided by using the freedom in renumbering the sites k — > k + ho, ko G Z. The 
consequence of this freedom is that, for any QW on Z, there are infinitely many 
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orders of the basis giving a CMV-shape transition matrix. Just as good as the initial 
order would be to take 

|fcot), l*o " H), |*o - It), \koi), \ko + It), |*o - 2|), \ko - 2f), |*o + H), • • • 

where ko is an arbitrary integer. The order of the basis given originally in (pQ) for 
QWs on Z can be understood as a folding of Z at site 0, so these other possibilities 
correspond to foldings at an arbitrary site ho. 

These new foldings lead to different CMV matrices, measures and OLP, any of 
them could be used to study a QW on Z. Since the presence of localization in a 
QW has to do with the Lebesgue decomposition of the measure, it is important to 
know how the measure changes with the renumbering of the sites. This is answered 
by the following result. 

Lemma 3.3. Given a QW on Z, the measures /j., fi corresponding to different 
foldings are related by 

dfi(z) = A(z)dfj,(z)A(z)\ 
where A is a 2 x 2-matrix polynomial. 

Proof. The transition matrices U, U related to different foldings are representations 
of the same transition operator with respect to basis which only differ in the order. 
Thus they are related by conjugation with a permutation matrix II, i.e., U = LTCTTI. 
On the other hand, the KMcG formula ensures that 

{(U+zl)(U-zl)- 1 ) j , k = [ l^Xj^dfi^XkitY, 

Jt t — z 

where (*)j, k stands for the (j, fc)-th 2 x 2-block. Thus, the Carathodory function of 
the measure fi is given by 

= / ^ Mt) = ((U + zl)(U - zl)- l ) , , 
J j t — z 

and similarly for the Carathodory function F of ft. 

Each subindex k stands for a pair of indices which we will denote by k s , s = +, — . 
If II transforms the indices + and 0_ into j r and k s respectively, then Ili j o + = 5ij r , 
n ii0 _ = $i,k s , and 

F{z) = (Ilt(17 + zl)(U - ^i)^ 1 n) 0i o 

f((U + zl)(U - zl)- 1 )^ ((17 + zl)(U - zl)- 1 )^ 
~ \{(U + zl)(U - zl)-i) ksdr {{U + zl)(U - zl)^) ks>ks ) ~ 

where and X^ stand for the upper and lower row of Xk respectively. This 
proves the proposition with 

A(z) = z l (J^i^y , f° r some I > 0, 

since X r - and X s k are 2-vector Laurent polynomials. □ 
We are interested in the following consequence of the lemma above. 
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Corollary 3.4. Two measures of the same QW on Z with respect to different fold- 
ings have the same mass points, and the support of their absolutely continuous and 
singular continuous parts coincide. 

Proof. Let /x, fi be such measures. Then, dfi = AdfiA^ for some matrix polynomial 
A, which implies that fl({z}) = whenever fj,({z}) = 0. Since there must be 

another polynomial A such that dfi = AdjxA , we conclude that the mass points 
of /x and /i coincide. The rest of the assertions follow similarly from Proposition 13.31 
and the invariance of the absolutely continuous and singular character of a matrix 
measure on T under the transformation dv — > AdvA^ for any matrix polynomial 
A. □ 

Unless we state specifically a different folding, the measure of a QW on Z means 
for us that one related to the folding at site given in (pQ). Nevertheless, the previous 
corollary ensures that we can refer to some characteristics of the measure without 
indicating any folding because they are common for all of them. 

Besides exploiting different foldings, the general results about localization for 
QWs on Z also require the use of Wiener's theorem on the unit circle (see [341 
Theorem 12.4.7]): for any scalar measure /i on T 

1 N f 

jv™o 9 N + 1 \Vn\ 2 = ^2\Ki z })\ 2 , Vk= I z k dfi(z). 

Thus, /i has no mass points if and only if linijv_ 5 . 00 2Af 1 +1 J2n=-N \^n\ 2 = 0, which is 
satisfied in particular when lim^oo /i n = 0. The complex numbers fi n are known as 
the moments of the measure /i. 

The following theorem is the main result of this section. It gives an interpretation 
of the single state localization in terms of the measure of the QW. 

Theorem 3.5. Let a QW on Z or Z + with transition matrix U and measure fi. 

(a) // /j, is absolutely continuous, no state exhibits localization. 

(b) If fi has a mass point, all the states \^S„g) exhibit localization except at most 
one state at each site k which must have a, (3 ^ 0. The existence of such a 
non localized state is mandatory when the singular part of the measure is a 
single mass point in Z + or a single pair of opposite mass points in Z. 

(c) The states \^f) which do not exhibit localization must satisfy 

^(z)fj,({z}) = 0, WzeT. (9) 

(d) // \i has no singular continuous part: 

(i) No state exhibits localization fi has no mass points -v^ 
<=> lim ^U n ^ = for all if), # e L 2 (Z). 

n— >oo 

(ii) does not exhibit localization <^ ([5j) <^ lim ipU 71 ^ = 0. 

n— ¥oo 

Proof. Statement (a) follows directly from Lemma 13.21 

1^) = l^cfia) does n °t exhibit localization if and only if ]xm n ^. 00 ipU n (x/j ( ^Qy = 
\im n ^ 00 ipU n (ipQ k ly = 0, which obviously implies that \im n ^ 00 ipU n ip^ = 0, i.e, 



lim / z n ^{z)dii{z)il)(zy = 0. 

n— >oo Jy 
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In other words, the n-th moment of the scalar measure dfi^ = t/»cZ/lit/»' converges 
to zero as n — > oo. Wiener's theorem ensures that has no mass points, which 
means that ip(z)fj,({z})tp(zy = for any z G T. Bearing in mind that [i({z}) is 
positive semidefmite, we get (c). 

According to (JSJ), given a QW on Z + , condition (EJ) becomes ift(z) = «I 2 t( 2 ) + 
(3X2k+i{z) = for any mass point z of fi. Since the OLP have no zeros on T, when 
fj, has a mass point this equation has a one-dimensional subspace of solutions (a, (3), 
a, (3 7^ 0, which represent the same quantum state. The existence of more than 
one mass point or a singular continuous part of the measure can give incompatible 
equations for a, (3, so the presence of states which do not exhibit localization is only 
ensured in the case of at most one mass point. This proves (b) for Z + . 

For a QW on Z, (JS} implies that the no localization condition at site k = 
reads as ij>(z)n({z}) = (aXj + f3X^)n({z}) = for any mass point z of /x, with 
the upper and lower rows of Xk- We know that Xo = 1, while Xi can be 
calculated from the first two 2 x 2-block equations of ([3]), 



On the other hand, the mass of any mass point z is the singular matrix given in 
Proposition 12.31 Combining these results we find that 



The coins are unitary, so | det Cq\ = 1. Besides, the assumption of the irreducibility 
for the QW implies that cfj ^ 0, so |c^| 2 = 1 — \ c fj\ 2 < 1- Since |^(^)| = 1 for 
any mass point z, the above equation becomes a + /3k(z) = with k(z) ^ 0. This 
equation is invariant under the reflection z —> —z due to Corollary 12.41 Therefore, if 
there is a single pair of opposite mass points, such an equation has a one- dimensional 
subspace of solutions (a, (3), a, (3 ^ 0, which represent the same quantum state. The 
presence at site k = of non localized states is ensured only when the singular part 
consists in at most a single pair of mass points, otherwise incompatibilities can 
appear between different equations for (a, (3). 

To generalize this results for k ^ is enough to use the folding at site k. Corollary 
I3.4l states that the mass points will not change when choosing this new folding. Thus, 
the previous discussion remains unchanged, but the conclusions are now about the 
states at the site k, which plays the role of the origin with this new folding. Therefore 
(b) is proved for Z too. 

Let us return to the general case of a QW on Z or Z + , and assume that the 
measure has no singular continuous part, which will be true no matter which folding 
we choose in Z, again due to Corollary 13.41 Then (a) and (b) are the only options. 
In case (a), Lemma 13.21 states that ]im n ^. 00 ipU n ^ = for all i^^ip. Also, as we 




giving 




(10) 



1>(z)p({z}) = & a + J^( c ^> + Z T]{z) det C ) = 0. 



16 



CGMV 



pointed out at the beginning of the proof, the condition lim^oo if;U n ^ = is not 
only a consequence of the fact that | ty) does not exhibit localization, but also yields 
([9]). On the other hand, ([9]) gives J T ip(z)d/j, s (z)^j}(zy = for any i/> because the 
mass points constitute all the singular part of the measure, thus implying that \^) 
does not exhibit localization. This finishes the proof of (d). □ 

The previous theorem indicates that the more mass points the measure exhibits, 
the less possibilities for non localized states because, apart from the conditions 
associated with the singular continuous part of the measure, there are as many no 
localization equations Q as mass points, which makes it more difficult to have non 
localized states as the number of mass points increases. 

The proof of Theorem I3.5I shows that, in Z, the no localization equation fl9]) is 
invariant under the reflection z —> —z, hence, only one of such equations must be 
taken into account for each pair of opposite mass points. 

There are situations in which it is known that the singular part of the measure 
is not purely continuous (or even no singular continuous part appears). Theorem 
I3.5I provides in such a case a localization dichotomy: either no mass points and 
no localized states exist, or there are mass points and "almost" any state (at most 
one exception per site) exhibits localization. When this dichotomy works we can 
talk about QWs with or without localization because then localization becomes an 
"almost" global property. 

Moreover, Lemma I3.2I shows that in the absence of a singular continuous part of 
the measure the asymptotic return probability can be computed exactly through 



zgT 



Ztz" 



where the sums are in fact over the mass points z of /x. 



3.1. Periodic QWs with finite defects. Among the QWs where the localization 
dichotomy works are those with periodic coins, with or without a finite number of 
defects. 

Proposition 3.6. // the coins Ck of a QW on 7* or Z + satisfy Ck+ P = Ck, pGN, 
for all but a finite number of sites k, the corresponding measure has no singular 
continuous part and thus the localization dichotomy holds. 

Proof. Consider first the case of strictly periodic coins on Z + with period p, i.e., 
Ck+p = Ck for all fc£Z + . The related measure \x has Verblunsky coefficients 

-(o) -(fc) -i( CT (°>+-+<T( fe - 1 )) n l \ 1 

a = c\{, a 2k = 4i e 1 j , a 2k -i = 0, k > 1, 

where = + ■ The fact that and have period p ensures that the 
new Verblunsky coefficients 

a k = a k e^\ = ±(a^ + ... + a (p-D) (i 2) 

2p 

have period 2p. 

Let fi and / be the measure and Schur function associated with the Schur pa- 
rameters dfc. As a consequence of the periodicity of &k, the Schur iterate f 2p has 
the same Schur parameters &k as /, hence f 2p = f ■ Bearing in mind that any 
step of the Schur algorithm (jSJ) is a rational transformation, the relation f 2p = f 
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can be written as a polynomial equation for f(z) with polynomial coefficients in z. 
Therefore f(z), and thus zf(z), are algebraic functions of z, which implies that the 
equation zf(z) = 1 has a finite number of roots. This means that the singular part 
of fi has a finite support, so it can not have a continuous part. 

Relation (1T21 between and 6>k implies that the corresponding measures /i, fx 
are connected by a rotation (see page 473]), dfi(z) = dfi(e~^z), thus /i has no 
singular continuous part neither. Besides, from the link between the measure fi and 
its Schur function / we find that f(z) = e %<d f (e _M? z) , thus / is algebraic too. 

Now suppose that we modify a finite number of coins Ck, so that Ck+ P = Ck 
only holds for k > k . Then, the sequence (ak)k>k is periodic with period 2p, 
and the corresponding Schur function, which is fk , must be algebraic. The Schur 
function / is obtained from ff. by ko steps of the inverse Schur algorithm (jSJ), each of 
them preserving the algebraic character. Hence, / and / are algebraic too, and the 
measures fi and \x have no singular continuous part, just as in the strictly periodic 
case. 

With regard to QWs on Z, the periodicity of the coins Ck for any k G Z with 
\k\ > ko implies again the periodicity of the Schur parameters &k given in (|T2|) 
for the same range of indices. Therefore, the previous arguments show that the 
Schur functions /+, /_ associated respectively with the Schur parameters a£ = ctfc, 
= — a-k-2 are algebraic. Since the singular part of the matrix measure /j. of the 
QW is supported on the roots of z 2 f + (z)f-(z) = 1, the result follows from the fact 
that z 2 f + (z)f-(z) is algebraic. □ 

In the case of QWs on Z with strictly periodic coins, stronger results can be 
achieved. 

Proposition 3.7. Any QW on Z with strictly periodic coins is free of localized 
states. 

Proof. If a QW on Z has strictly periodic coins, the full sequence (d fc ) fcgZ appearing 
in the proof of the previous proposition is periodic too. The matrix measure fi of 
the QW is a rotation of the measure fi with Verblunsky coefficients 

\&k J 

The block CMV matrix C with Verblunsky coefficients ctk is obtained by folding 
a two-sided CMV matrix C with scalar Verblunsky coefficients &k (see [5]). Like any 
two-sided CMV matrix with periodic Verblunsky coefficients, C has an absolutely 
continuous spectrum (see [7]), and the same holds for C because it is related to 
C by a simple reordering of the basis. This means that the scalar measure fi^ 
defined by tjjC ^ = J T z n dp, 1 p(z), n e Z, is absolutely continuous for all ip. Then, 

lim n _ 5 . 00 ip C ip^ = for any ip and Theorem I3.51 d.ii implies that no state exhibits 
localization. □ 

3.2. Quasi-deterministic QWs. QWs with diagonal coins Ck for any k are de- 
terministic because the one-step transitions 

Z |jfet) |fc+lt) \kl)^\k-li) 

Z+ > \2i) -»■ |H) |0|) |0|) 1 1 T> -»■ |2|) -»• ■ ■ • 
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take place with probability one. These QWs exhibit no localization, even when a 
finite number of defects appear, regardless of the number and details of the defective 
coins. 

The presence of a finite number of defects means that C k is diagonal for all but a 
finite number of sites k. In such a case, the related measure has null Verblunsky 
coefficients ct k except for a finite number of indices k. Hence, ct k = for k > ko 
and is a Bernstein-Szego measure which can be expressed using the OLP as (see 
[8] for the general matrix case) 

d^) = [X ka {e^X k0 {^)]-^. 

Since \x is absolutely continuous, no state exhibits localization. 

Concerning the possibility of having a singular continuous part in the measure, it 
is known that sparse sequences of Verblunsky coefficients on Z + can give a measure 
which is exclusively singular continuous (see [11] and [3H Section 12.5]). This shows 
that such a pathological situation can appear surprisingly close to the deterministic 
case corresponding to diagonal coins. 

Another source of singular continuous measures are those measures supported on 
a Cantor type set (see [10] and [331 Section 2.12]) or those given by appropriate 
infinite Riesz products (see [32] and [331 Section 2.11]). It would be interesting to 
search for QWs corresponding to these kinds of measures, as well as to study the 
localization properties of such rather pathological situations. This could shed light 
on the general picture for the localization properties of QWs with a measure whose 
singular part is purely continuous. 

The models which we will analyze in detail are the QWs with a constant coin up to 
one defect at the origin. They are a special case of periodic QWs with finite defects, 
so the localization dichotomy works for them. We will make an exhaustive analysis 
of localization in these examples, both on Z and Z + , to illustrate the effectiveness 
the CGMV method beyond the case of a constant coin, and to understand how the 
single state localization depends on the parameters of the models. 

4. QWS WITH ONE DEFECT 

We will consider a general QW with coins C k which are constant except for the 
site k = 0, i.e., 

C k = C= ( Cu ° 12 ) , k^O, C = D=( d , n d , 12 ) . (13) 

\C21 C22 J ^"21 "22 / 

We will refer to this as a QW with one defect at the origin. To place the defect at 
the origin is simply a convention for the numbering of the sites in Z, but it is real 
restriction in Z + . 

Remember that we only need to consider irreducible QWs, which means that we 
can assume without loss of generality that Cjj, djj 7^ 0. We will use the notation 

a a 

a = a 1 + a 2 , r = n + r 2 , & = -, 

where e lcrj and e lTj are the phases of Cjj and djj respectively. 

The case of a diagonal coin C is somewhat special. We know that it leads to an 
absolutely continuous Bernstein-Szego measure which, therefore, yields a QW with 
no localization. Thus, for the general discussion we will suppose that c 2 i 7^ 0. 
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4.1. QWs with one defect on the non negative integers. Let us consider 
first the coins f[T3"j) on Z + . According to the results described in Section [21 the order 
indicated in (CQ) gives a transition matrix U = ACA^ , A = diag(l, Ai, A2, • • • ), with 

A 2fc _ x = e i(™+(fc-i)«) j A 2fc = e -*Cn+(fc-i)<xi) ? fc > Xj 
and C = C(o!fe) a CMV matrix with Verblunsky coefficients 

)d 2 i, i£k = 0, 

a2k = \c 21 e-^-^\ if*>0, a2fe+1 = °' k ~°- 

The Verblunsky coefficients can be written as = ake~ z ^ k+1 ^ with 

(d fc ) = (6,0,a,0,a,0,a,0,...), a = c 2l e i(|<J ~ r) , b = d 2 ie^. (14) 

This means that the measure //, the Carathodory function F and the OLP of 
the QW are related to those ones of C = C(ak) by (see [5j page 473]) 

F(z)=F(e- w z), X fc (z) = A fc X fc (e- w z), (15) 




A 2fe _i = A 2fc _ ie -^ = e ^^+r 2 -, 2)) ^ ^ ^ 
A 2fc = A 2fc e* M = e^- 



with an obvious notation for the elements corresponding to C. In other words, the 
OLP of the QW are, up to a change of phases, a rotation by an angle 1? of those 
corresponding to a CMV matrix with Verblunsky coefficients (b, 0, a, 0, a, 0, a, 0, ... ). 

The related Schur function / = f a ^ has Schur parameters (6, 0, a, 0, a, 0, a, 0, ... ). 
Its second Schur iterate f 2 is the Schur function f a whose Schur parameters are 
(a, 0, a, 0, a, 0, ... ), so from (jSJ) we find the relations 

, / x 1 fa,b( Z )~ b f ,\ Z 2 f a (z)+b 

Uz) = W*> = YTWuT) (16) 

In particular, setting b = a, f a ,b becomes f a . This leads to the quadratic equation 
az 2 f a (z) 2 + (1 — z 2 )f a (z) — a = for f a which yields the expression 



/.(*) = ^ 1 ^ K ^ , ^ = ^ ~ !) 2 + (17) 



Since f a is analytic in D, the choice for the square root must result in the branch 
such that a/ A a (z) z ~*°> 1. Such a choice implies that the boundary values of f a on 



the unit circle ar 



Ra{0) 



-id 

f a (e w ) = —(R a (8) + 1 sin 9), 
a 

sgn(cos#)\/|a| 2 — sin 2 9, if |sin#| < \a\, 
■isgn(sm9)^/ sin 2 9 — \a\ 2 , if |sin#| > \a\. 



Thus, \fa(e i8 )\ = 1 if | sin 6> | < |a| and \f a (e ie )\ < 1 if | sin6>| > \a\. This also 
holds for / = / 0> 6 because any step of the Schur algorithm preserves the relations 
< 1 and |/(2)| = 1 at any point z G T. Therefore, according to @, the 
weight w(9) of the measure \x lives on |sin#| > \a\, which defines two arcs which 
are symmetric with respect to the real axis. The singular part is supported on the 



2 See [5] Appendix] for a discussion about the boundary values of y/A a on T. 
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finite number of roots z G T of zf(z) = 1, so it can have only mass points. When 
these mass points are present, they must lie on any of the two complementary arcs 
given by 

T a = {e ie : |sin0| < \a\}, 

which are symmetric with respect to the imaginary axis. This is because the equality 
zf(z) = 1 implies |/(^)| = 1 for any z6t The consequences of these conclusions 
for the measure of the QW are obvious because \i is obtained simply rotating fi by 
an angle 

Different coins f|T3|) giving the same pair a, b G D have measures which only differ 
in a rotation. Concerning the relative values of a and b, when the defect disappears, 
i.e., D = C, we get b = a = C2\e 1 ^ . Nevertheless, the defect not only changes the 
first Schur parameter from a to b, but it affects also the value of a which acquires an 
extra phase e^ a ~ T \ Due to this, the equality b = a can happen even with D ^ C, 
indeed it is equivalent to d 2 \ = C2ie t<yT ~ a \ The restriction c 2 i ^ that excludes the 
special case of a diagonal coin C means that we are considering a ^ 0. 

4.2. QWs with one defect on the integers. Assume now that we have the coins 
(TT5j) in Z. From the general results of Section [2] we find that the order indicated in 
dl]) gives a transition matrix U = ACA 1 * , A = diag(l, A 1; A 2 , . . . ), where 



ikai q \ /g-i(n+(fc-l)oi) q 



*2fc-l ~ \y e i(T2 + (fe-l)<7 2 ) J' ^2fc - ^ Q e -ika 2 J , > 1. 

and C = C(ajfc) is the CMV matrix with Verblunsky coefficients 

/ -«_ 2 fc-2\ n 7 \ n 

OL2k = I n , «2fc+l = 0, k>0, 



U 2 k 

{d 2 i, if k — 0, 

C 2ie -i(r+(fe-l)a) ) if Jfc > 0, 
c 2 ie~ ika , if Jfc < 0, 

As in the case of Z+, a rotation plays a useful role in our analysis. Defining, 



•I I i ^ ~ r ■ C 21 ^— — 7 • C 21 i(--a') / 1n \ 

z j c 2 i | e 2 , 6 = 2- re 2 d 2 l, W=Z-j re H 2 , (19) 

l C 2l| l c 2l| 

we can rewrite CKfc = e _ ^ fc+1 ^o;fc, where 

(d fc ) = (/3, 0, a, 0, a, 0, a, 0, ...), a = ^ ^ a J , /3 = Q? & ^ a V 

As a consequence, the measure ^, the Carathodory function F and the OLP 
of the QW are given by 

= dii{e- u z), F(z) = F(e-^z), X k (z) = X k X k (e~ l{> z) , (20) 



1, 



e 



V2/c — ^2/c 



fc > 1, 



(TJ-CT2 
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where elements with a hat are related to C — C(ctk). Therefore, just as in the 
case of Z + , up to phases, the OLP of the QW are obtained rotating by an an- 
gle $ the OLP going along with a one defect sequence of Verblunsky coefficients 
((3, 0, a, 0, a, 0, a, 0, . . . ). We should remark that the defect is only at the (2,l)-th 
coefficient of j3. 

The corresponding Schur function / has as Schur parameters the antidiagonal 
sequence (/3, 0, a, 0, at, 0, a, 0, . . . ). Applying Propositions 12.31 and 12.21 we conclude 
that 

where f a and f a ^ are the scalar Schur functions introduced in the previous sub- 
section, although here a and b bear a different relation to the coefficients of the 
coins. 

The measure /j, of the QW is obtained as a simple rotation of ji, so we only need 
to discuss this last one. We know that |/ a | = \f a ,b\ = 1 in the two closed arcs Y a 
and | f a |, \f a ,b\ < 1 in the two open arcs T\ T a . Therefore, the same result holds for 
ll/H = max{|/ a |, |/ a ,b|}- Using (j3J), we find that, for a.e. 9, the weight w{9) of fi is 
singular if and only if e ld G T a . Furthermore, (HI) also yields for a.e. 9 

= (1 - e- id f)-\l - /(e*)t/(e«))(l - e*/)- 1 . 



The equality 



ti_A-|/a, 6 | 2 



l-|/a| 2 

shows that detw(9) = implies \f a ,b{e. l6 )\ = 1 or \f a ,b{^ e )\ = 1- Since these two 
conditions hold simultaneously, detu>(#) = necessarily gives f(e %e )'f(e ie ) = 1 
and thus w(9) = 0. We conclude that w is zero in r a and non singular in T \ r a . 

As for the singular part of fi, it is supported on a finite number of points, the roots 
z G T of z 2 f a (z) f a ,b(z) = 1, and any of these roots must satisfy \f a (z)f a ,b( z )\ = 1- 
Hence, the singular part only can have mass points located at T a . 

In contrast to the case of Z + , three parameters a, b G D, oj G T characterize now 
the coins ( Tl3l) with the same measure up to rotations. On the other hand, just 
as in the case of Z + , the equality b = a does not hold only for D = C because, 
remarkably, it is equivalent to the same condition o^i = C2\e l ^ T ~^ appearing for the 
non negative integers. Also, the consequences of the defect are not only encoded in 
b, but the imaginary value a = i\c2i\ for a constant coin C acquires with the defect 
an extra phase e tS ^~ which is the square root of the similar extra phase for the case 
of Z + . As in Z + , we only need to consider a^O because we know that a = yields 
no localization. 

5. Localization: one defect on Z 

The previous discussions indicate that the study of localization in QWs leads to 
the analysis of the mass points (in general, the singular part) of the corresponding 
measure. This analysis is worth doing specially for QWs where the localization 
dichotomy works, such as periodic QWs with a finite number of perturbations. The 
QWs with one defect are just the simplest examples of this case. As we pointed 
out, in such situations we will talk about QWs with or without localization because 
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then localization can be viewed as a global property: it holds for no state or for 
almost any state. 

Surprisingly, for QWs with one defect, the mass points analysis is simpler for Z 
than for Z+, among other reasons, due to the symmetry of the mass points with 
respect to the origin, which is lost for Z+. Hence, we will study first localization in 
a QW on Z with coins (TT3j) . 

Concerning previous related results, N. Konno has proved in [26] that the state 
-^|0 t) + ^75 1 D exhibits localization in Z for the perturbation of the constant 
Hadamard coin 

given by 

C = H ' D = 7i i?-* -i)- < 21 > 

whenever e 1 ^ 7^ 1. We will recover this result as a particular case of our analysis but, 
assuming it for the moment, observe that we have a stronger result: the dichotomy 
implies that all the states should exhibit localization up to, at most, one state 
per site. Indeed, we will see that in this model any state cn|0t) + /?|0i) exhibits 
localization. 

The aim of this section is to perform a systematic study of localization for any 
QW with one defect on Z, which can reveal in the simplest examples the coin 
dependence of localization properties. 

Therefore, our purpose is to determine which coins ffTSl give in Z a measure with 
mass points. Subsection 14.21 shows that these models fall into groups with a common 
measure up to rotations, each such a group characterized by the three parameters 
a, b G D, co G T given in ( fl9l) . A canonical representative of the measures in a given 
group is that one fi — fx^ b associated with the common CMV matrix C = C(a^) of 
the group given in Subsection 14.21 whose weight and mass points live in T \ T a and 
T a respectively. 

Coins (I13p with the same values of these parameters have the same mass points 
up to rotations and, therefore, the corresponding QWs have the same localization 
character. For instance, in Z, the constant Hadamard coin has the same values 
a = b = -4=, co = 1 as its perturbation 

1 f e i<t> 1 
C = H, D = — ( 

which proves in a very simple way that no state in such a QW exhibits localization 
(see [25] for a different proof in the particular case of the state ^|0f) + 75IOI)). 

5.1. Mass points of fi^b- Bearing in mind that localization in a QW with one 
defect on Z only depends on the associated parameters a,b,co, we can restrict our 
attention to the canonical representative /j% b . Proposition 12.31 states that the cor- 
responding mass points are the roots z G T of g a ,b{z) = z 2 f a (z) f a ,b(z) = 1 such 
that 

m a , b (z) = lim - — ± 0. (22) 

These conditions do not depend on co, so the mass points, and thus the localization 
behaviour, only depend on a, b but not on co. 
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The choice of the square root y/A^ makes f a analytic in T except at the branch 
points, i.e., the solutions of A a = 0, which are the four boundary points dT a of the 
two arcs T a , 

dT a = {±z a ,±z a }, z a = p a + i\a\, p a = \J\ - \a\ 2 . 

The relation (1TB|) between f a and f a b shows that f a ^ is analytic in T \ dY a too, and 
so the same is true for g ab . In consequence, we find from Proposition 12.31 that the 
measure fi^ b has a mass point at any root z G T \ dT a of g a ,b(z) = 1- Indeed, we 
know that these roots must be in the interior T° of T a because there is no root in 

T\r . 

Therefore, the roots z G T of g a ,b{z) = 1 can lie on r°, and then they are mass 
points for sure, or they can be on dT a , in which case we should check ([22]) to decide 
if they are mass points or not. 



5.1.1. Mass points on dT a . We will prove that, although the points of <9r a can be 
roots of g a ,b( z ) = 1) they are never mass points of /z£ 6 because condition (1221) is not 
satisfied. We will only consider z a , the analysis for the remaining three points is 
similar. 

First, let us find the values of b which make z a a root of g a ,b{z) — 1- We know 
that |/o(2 )| = 1 because z a G T a , so from (1161) we obtain 

g a , b (z a ) = z2 M^ + \ 

and g a ,b(z a ) = 1 becomes equivalent to Jm(z^f a (z a ) + b) = 0. On the other hand, 
the expression ([17]) for f a gives z1f a (z a ) = i&z a , so that z a is a root of g a {z) = 1 if 
and only if Ivab = — lm(iSvZ a ). 

Now, assume that b satisfies the condition Imb = — Irn(z-A-z a ), and let us compute 
linvj! m(rz a ). The first order Taylor expansion of A a (rz a ) at r = 1 gives 

A a (rz a ) = K^l - r) + 0((1 - r) 2 ), K x ± 0, 

which, using ([IT]) , yields 

(rza) 2 fa(rza) = i-^z a + K 2 JT^ + 0(1 - r), K 2 ^ 0. (23) 
\a\ 

Inserting this into the relation 

M , _ (^(j + l)(z 2 f a (z) - 1) + ^(Im6)z 2 / a (z) 

1 + bz 2 f a (z) 

obtained from ([IE]) , leads to 

g a ,b{rza) ~ 1 = KRe{i±z a )y/T=? + 0(1 - r), AT ^ 0. 

This proves that lim^ m(rz a ) = when Re(iAz a ) ^ 0. The equality Re(i A z a ) = 
would imply |Imfe| = |Re(z-A-z a )| = 1, which is not possible, thus we conclude that 
there is no mass point at z a , even if it is a root of g a ,b(z) = 1. 
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FIGURE 1. The transformation z i— > ( maps both and T~ one 
to one onto E a . 5(a) is the open set limited by the arc S a and the 
straight line passing through and m S re y color in the figure. 



5.1.2. Mass points on T®- At this point we know that the mass points of fx^ b are the 
roots of g a ,b{ z ) = 1 in We can restrict our analysis to the right arc r+ = {e 10 G 
r° : cos 9 > 0} of T° because the mass points appear in pairs ±z, one belonging to 
r+ and the opposite one lying on the left arc T~ = {e ld G T° a : cos# < 0}. 
The study of the roots in r+ is simplified under the change of variables 

C = ((z) = -z 2 f a (z), 

which maps r+ one to one onto the arc (see figure ED) 

£ a = {|^:cost< \a\} : 

and dT+ = {z a ,z a } onto <9£ a = {C,C}, Ct = ± ipa)- These mapping 

properties can be inferred from the expression 

46 



((e w ) = - e — ( \j\a\ 2 - | sin 2 6>| + *sin# ) , e ie G T+ 



a 



obtained from ( TT5|) . which shows that ((z a ) = ( a , C(^a) = Co*" an d the argument of 
({e ie ) is increasing in 9 for e* 61 G T+. The inverse mapping is 

2 = z C) = — =7,- 
|1 - aC| 

The arc E a can be alternatively described as 

S a = {( G T : Re(aC) < |a| 2 } = {C G T : \a - || < ±}, 

a result that will be of interest later on. 

Bearing in mind that \f a \ = 1 in r+, we find that 

so the translation of the equation for z to the new variable ( is 

g a>b (z) = 1, z G T+ ^ Im& = ImC, C e E a . 
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Given b G D, the solutions ( G T of the equation Imb = Im( are 

C±(6) = ±Vl~ lm 2 b + ilmb. 

The values of a which are compatible with (±(b) are given respectively by any of 
the equivalent conditions 

C±(6)GS a & Re(aC±(6)) < |a| 2 & \a - §C±(&)| > \- (M±) 

Therefore, the measure has mass points if and only if at least one of the 
conditions M + , M_ is satisfied. For each of the conditions M + , M_ which is 
satisfied, there is a pair of mass points at ±z + (a,b), ±Z-(a,b) respectively, where 

z±(a,b)= )~^f t l GT+ -z ± (a,b)eT-. (24) 
|1 - a(±(b)\ 

Hence, we have the following possibilities: 

(M°) If none of M-t are satisfied, has no mass point. 

(M 2 ^) If M + is satisfied but M_ is not, /x" 6 has 2 mass points: 

z + (a,b) G T+ and — z + (a, b) G Y~. 
(M 2 ) If M„ is satisfied but M + is not, /x^ 6 has 2 mass points: 

Z-(a, b) G T+ and — z_(a, b) G T~. 
(M 4 ) If M-i- are both satisfied, fi^b nas 4 mass points: 

z±(a, b) G T+ and — z±(a, b) G r~. 

The case of 2 mass points is characterized by M 2 = (M+ or M 2 ), while M 
YI . or M_) is the condition for the existence of mass points. 

5.2. Localization pictures on Z: dependence on a and b. The localization 
dichotomy for one defect on Z does not depend on u, but only on a, b. Hence, we 
can discuss two kind of problems: Given a, which values of b yield localization? 
Given b, which values of a yield localization? 

5.2.1. From b to a. The last way of expressing M± in (M±) above means that a 
lies outside the closed disk Vf of center C±{b)/2 and radius 1/2. Therefore, the 
different cases can be stated as (see figure [2]): 

(M°) a G V~l H /i£ 6 has no mass point. 

(M+ ) a G T>^ \ T>^ -vv- /z£ 6 has 2 mass points ±z + (a, b). 

(M 2 ) aeD+\P^ j*£ 6 has 2 mass points =bz_(a, 6). 

(M 4 ) a^D+UD^ /*£ 6 has 4 mass points ±z+(a, 6), ±z_(a, 6). 
We conclude that, given a value of 6, a QW with one defect on Z exhibits localization 
if and only if a <£ V£ n XV . 

5.2.2. From a to b. Now we look at the first condition in (M±) above. To make it 
more explicit let us decompose S a = £+ U £~ into its right and left parts = 
{( G S a : ReC > 0}, £~ = {( 6 S a : Re£ < 0}. Then, a choice of a fixes E„ , and so 
the interval ImE„ where Imb must lie to fulfill M±. This means that localization 
for QWs with one defect on Z only depends on a and Imb, but not on Reo. 

Taking into account that the angular amplitude of E a is bigger than ir (remember 
that we are considering a ^ 0), we find three possibilities for Im£ a (see figure |3]): 

(A) (ReO(ReC) > & (-1,1) c Im£ a . 
(B+) (ReO(ReC) < 0, Ima > & ImS a = [-1, r+), r+ < 1. 
(B_) (ReO(ReC) < 0, Ima < ImS a = (r_, 1], r_ > -1. 
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Figure 2. Localization for one defect on Z (from b to a). In 

green color the values of a giving localization for the choice of b in 
red. They only depend on Imb. In light green the values of a with 2 
mass points and in dark green those with 4 mass points. 

The value (Re£+)(Re£~) = ^ ^Jp a turns these three cases into the following local- 
ization criteria: 

(A) Localization Mb <^> |Ima| < |a| 2 \a — || > \ or \a + || > \. 
(B + ) Localization for Im& < r + £ (0, 1) -vv- Ima > |a| 2 |a — || < |. 
(B_) Localization for Imfe > r_ £ (—1,0) Ima < — |a| 2 |a + || < \- 

Roughly speaking, the values of a split into three regions delimited by two circles 
with radius 1/2 centered at ±z'/2. Outside these circles localization holds for any 
defect. Inside the upper or lower circle there is respectively an upper and lower 
bound for Imb which delimits the defects giving localization (see figure E]). 

Notice that, for any a e D, localization holds at least for b lying on the open set 
S(a) limited by the arc S a and the straight line joining an d (a ( see figured]), 
that is, 

S(a) = {r-^e u : cost < \a\, r < 1}. (25) 

Indeed, S(a) yields exactly the values of b giving localization when a is imaginary 
because in that case Im£+ = Im(^. Thus, among the values of a with the same 
modulus, the biggest region of values of b without localization holds for Rea = 0. 
Since ImS+ = ImS~ for an imaginary value of a, it also ensures 4 mass points in 
case of localization. 

For a fixed a, the bounds r± = r±(a) are r + (a) = maxlmj^, (~} and r_(a) = 
minIm{C+,C}, i-e., 

r±(a) = Ima ± -^-|Rea|. 

\a\ 

These bounds also permit to distinguish between the values of b giving 2 or 4 mass 
points, once a is chosen. There are 4 mass points when Imb £ ImS+ D ImE~, 
and only 2 mass points if Im& £ ImS+ \ ImS^ or Im6 £ ImS~ \ ImS+. Looking 
separately at the three previous possibilities we find that 2 mass points appear when 
b lies on a band limited by two horizontal lines passing through and (~ . Hence, 
the situation in the three cases above can be more precisely described as follows 
(see figure [3]) : 
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Figure 3. Localization for one defect on Z (from a to b). The 

upper figure shows in dark green the values of a giving localization 
for any b. The upper (lower) circle in light green are the values of 
a such that localization fails for b lying on an upper (lower) band 
Imb > r + (Imb < r_). The lower figures represent in red color the 
values of b giving localization for each of the three values of a shown 
in the upper figure. They are characterized by Imb e ImE a . A pair 
of mass points appears for each of the conditions Imb e which is 
satisfied. Therefore, the dark red covers the values of b with 4 mass 
points, while the light red covers those with 2 mass points. 



Imb < r_(a), 4 mass points, 

(A) ^ r_(a) < Imb < r + (a), 2 mass points, 

r + (a) < Im&, 4 mass points. 

Imb < r_(a), 4 mass points, 

(B + ) ^ r_(a) < Imb < r + (a), 2 mass points, 

r + (a) < Imb, no mass points. 

Im6 < r_(a), no mass points, 

(B_) ^ r„(a) < Imb < r + (a), 2 mass points, 

r + (a) < Im6, 4 mass points. 



6. Asymptotic return probabilities: one defect on Z 

To compute the asymptotics of p£ l(n) as in (II ip we need, not only the mass points 
(known from the previous results), but also their masses and the OLP related to the 
site k. Let us see how to make the computations with the canonical representative 
dfi = d\x UJ ah instead of the actual measure d/j,(z) = dfi(e~ l ^z) of the QW. 



28 CGMV 

* ha 

Introducing f )20|) in ([8]) we obtain the relation il> a p(z) = ip & p(e~ l ®z) between the 
corresponding functions for the state l^^l), where 

a = \^a, (3 = Ag+i/5, if k = j, j > 

a = \ { 2j +1 a, = A^/3, if k = -j - 1, 

and A fc = di&g(X^\ A^). In particular, ^>i*o(z) = /^^^(e - "^) and tp^J(z) = 



K j ^ifli 6 % ^ z ) with Kj E T. Hence, fTTTj) can be written as 



2 



^z>l^( 2 )A({4)V'l fc dw t + 



2 



(26) 



For convenience, while performing the calculations we will omit the hat on d, (3 so 
that at the end of the computations we should make the substitution a, (3 — > a, (3. 
We also remember that p^ f3 (2n — 1) = 0, thus we only must consider p^ j3 (2n). 

6.1. Masses of fi = /x^ 6 . There are 4 possible mass points: zhz + (a,b), ihz_(a,b). 
We only need to calculate the mass of the two points z±(a,b) given in (1241) be- 
cause the mass of the opposite points follow from Corollary 12.41 We will make the 
calculations for a general point of the form 

zo = Co G E a . (27) 

|1 - aCol 

The mass of z±(a, b) is obtained setting ( = (±(b). 
Proposition 12.31 states that 

M {*o} = m(zo) -t-t 1 J ' 1 / / / \ 28 

\V{z ) 1 y m(z ) = l/z g ajb {zo). 

In this case 

rj(z ) = -ujz^Cq = -uj-r^ r, (29) 

ICo - a\ 

because we know that inverting (1271) yields Co — ~ z ofa( z o)- 

For the calculation of g' ab (zo), first perform the change ((z) = —z 2 f a {z) in g ai b(z), 

9a,b — =T, 

1 - b( 

so that 

9'a, b _ C L PU 

9a, b CV (l-6C)(C-6) 

and 

, , \ C(*b) A , > \ C / (^o)l-Re(R ) 



Co V IC0-&IV Co ICo — 

It only remains to compute C'( z o)- From ( fT7l) we obtain 

C(z) = 2z C -^, 

hence 

C'(zq) 1-aCn 1 1 — aCo I ,„ >x 
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1 |Co-&| 2 |a| 2 -Re(aCo) 1 1 " 
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2|Co-a| 2 l-Re(6Co 



(31) 
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6.2. Asymptotics of p^g(n) on ^- The asymptotic return probability to the origin 



(0) 



involves ip^'a = (a,0)X + (0,/3)Xi. We know that X = 1, while Xi can be 
obtained specializing the general expression ( TTOl) for the coins 

Pb -uP 



Pa Wffl 

-ua p a 



K ub p b 

related to the CMV matrix C of fi. This gives 

Z^/pa -Ua/p a 



Pb 



X 1 (z) 



-ub/p b z 1 /p b 



We finally find that 



^(z) = a(l,0) + z- 1 ). 

Pb 



6.2.1. The case of 2 mass points. Assume the case M 2 = (M^ or M 2 ), so that 
there are exactly 2 mass points: ±z + (a,b) for M+, and ±£_(a, b) for M 2 ,. For 
convenience, let us write z Q = z±(a, b) and Co = C( 2 o) for M±. Then, from ( )28|) and 
(El we find that 



V'l ^o)A({2o})V'^o(2o) t = m(zo) la - (3u 

pb 



(0) 



Pb 



0- C0 + b 
a — puj 

Pb 



where m(zn) is given in (l3~lj) . 

The result for — zq is the same because m(z) is even and rj(z) is odd, thus, ac- 
cording to (12"B1) . 

2 



lim pi°l(2n) = (2m(*o)) a ( 1 + 

Setting z = z±(a,b), then ( = C±(&) an d 

\Co+W 



2 A ICo + feh 



ct — /3w 



1 + 

which gives 
with 



ill 



IC0-&I 



1 + 



2V1 - Im 2 6 



^6 



V 7 ! -Im 2 6^Re6' 



lim^(2n)=^(a,&,u;) for M 2 , 



P 2 a 



1 



IC0-&I 2 

2 x 2 



a — /3cj 



Pb 



(32) 



lC±(6)-o| 



IT 



|d| 2 - |/3| 2 )Re6 + 2p fe Re(wd/3) 
V 7 ! - Im 2 6 
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Here we have used that \a\ 2 + |/3| 2 = 1. 

We see that the asymptotic return probability Ja,b,u) to the origin depends 
on the coefficients a and (3 of the state, as well as on the parameters a, b, u associated 
with the QW. Indeed, there is a state «|0t) + /?|(H) which exhibits no localization, 
given by 

Pb r„„ ™2 



a uj- 



for M 



Reb± VI - Inrfe 



(33) 



6.2.2. The case of 4 mass points. In the case M 4 there are 4 mass points: ±z + , 
where z + = z + (a,b) is related to ( + = C+(&)> an d i- 2 -) where z_ = z_(a,b) is 
related to £- = (-(&)■ Therefore 

2 



(0) (2n) 



2m(z) I a 

S=Z± ^ 



2m(z] ^ 



pu 

Pb 



. 2n 



z=z± 



a 



. 2n 



Pb \ Pb 

The cross terms of both summands cancel each other because 



— I, 



lim 3 l(2n) = p+Ja, b, u) + p a Ja, 6, w) for M 4 , 



hence, 



with p a p(a, 6, u) given in (|32|) . 

In other words, the 4 mass points ±£+, ±z_ contribute to the asymptotic return 
probability to the origin simply by adding the contributions that they should have 
if considered as two independent cases with 2 mass points. As a consequence, the 
existence of 4 mass points ensures that all the states at the origin exhibit localization 
because the two conditions in f l3"B"|) are incompatible for (a,f3) ^ (0,0). 

Particularly simple is the case of an imaginary value of a which, according to 
(fT!?]) . corresponds to a defect such that e %T = e l ° ". Then, localization appears if and 
only if Ima > 0, Imb or Ima < 0, Imb, and in such a case there exist always 4 mass 
points. The simplicity of the asymptotic return probability to the origin comes from 
the fact that, for an imaginary a, we have \( + — a\ = |£_ — a\ and thus 

/ 2 \ 2 /OT/T T7N\2 

,(0) fn„\ _ / 1 Pa 



lim Pa' 8 (2n) 



2Ima(Ima — Imb) 



(34) 



\(± — a\' 2 J \1 + lm z a — 2ImaImo, 

In this case the asymptotic return probability to the origin does not depend on the 
state. 

For instance, the model 02 ip gives a = -y=, b = ^=-, oj — 1, which exhibits 



V2' 



localization when Imb < Ima, i.e., e 1 ® ^ 1. The application of (|34|) to this model 
yields 



lim p^i{2n) 



2(1 — COS( 

3 — 2 cos ( 



which shows that the result obtained in [26] for the special case ct = ^75, /5 = ^75 is 
indeed true for any a, j3. 
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6.3. Maximum asymptotic return probabilities on Z. The previous results 
seem to indicate that the maximum values of 

Pa!{3 = lim PaU 2n ) 

should be reached for a close to the unit circle, which means that the non defec- 
tive coin C has an almost anti-diagonal shape. Let us analyze the behaviour of 
max aiJ gp^g when \a\ — > 1. 

Given a value of a with a fixed phase different from that of (±{b), the localization 
pictures in subsection 15.21 show that the measure has 4 mass points as far as \a\ 

is close enough to 1. Take a , A) such that f3 = iua , so |a | = \fto\ — Then, 
we find from ( 1321) that, in the case of 4 mass points, 

™fP a ,p > PaoJh ~ 2 [ l | C+ ( 6 ) _ a \*J + 2 V 1 |C_(6) - a| 2 J " 
Therefore, 

lim maxp^l = 1, a G T \ {(±{b)}. 

a— >ao a,p 

That is, if Ima ^ Imb and |a| is close enough to one, there exist qubits which 
asymptotically return to the origin with probability almost one. According to ( TT9]) . 
given a defect D, this holds for almost any coin C as long as its diagonal is close 
enough to zero. 

These results become stronger when a is imaginary. Then is independent 

of a,/3 and (1341) yields lim^-i-j = 1 for any state. In other words, if a is close 
enough to i or —i, all the qubits asymptotically return to the origin with probability 
almost one. Looking at ( Tl9l) we see that, given a defect D, this is the case of any coin 
C which is close enough to an anti-diagonal one provided that e ta is close enough 
to e iT . 



7. Localization: one defect on Z + 

We will study the localization for the coins f ll3p in Z + . As in the case of Z, this 
requires the analysis of the mass points of the corresponding measure. Subsection 14. II 
shows that these models fall again into groups with the same localization behaviour 
because any such a group has a unique measure up to rotations. Nevertheless, these 
groups are characterized now by only two parameters a,i G D given in fTl4l . The 
measure \x = of the CMV matrix C = C(dfc) introduced in Subsection 14. II serves 
as a canonical representative for the measures in a group. The corresponding weight 
and mass points are supported in T \ T a and T a respectively. 

7.1. Mass points of fi a ^. Concerning localization properties for one defect on 
Z + we can restrict ourselves to the measures n a ,b without loss of generality. The 
corresponding mass points are the roots z G T of h a ^(z) = zf a) b(z) = 1 such that 

= Iff i XArz) + °" (35) 

These roots must lie on T a and when they lie on L° condition f l3"5"j) is always satisfied 
because h a> b is analytic in L°. 
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Moreover, the points of dT a can be roots of h a< b(z) = 1 but never mass points of 
fi at b because condition (1551) is not satisfied on dT a . Consider for instance the point 
z a and assume that h a> b(z a ) = 1. Then, using (1231) we find that 

h{rz a ) - 1 = K'VT^r + 0(l -r), #V 0, 

which, according to (135]) . implies that /i a ,b({^a}) = 0. A similar proof works for the 
remaining points of dT a . 

Therefore, the mass points of [i a ,b are exactly the roots of h a ,b{z) = 1 in r°. To 
study these roots we will use the same chage of variables as in the case of Z. However, 
since the symmetry of the mass points with respect to the origin disappears in Z + , 
we must study independently the roots in the right and left arcs of r°. 

The transformation 

C = C(z) = -z 2 f a (z) 

maps both arcs onto S a (see figured]), with and inverse mapping given respec- 
tively by 

z = ±z{Q g r± .(C) 



Therefore, the equation for zG^ reads in terms of ( £ S a as 

h a , b (z) = l, zETi =1, (6E tt 

IC-6| 2 IC-a| C-a 

The last of the above equivalent conditions states that b lies on a straight line 
passing through £ in the direction given by i^C, — a or \/£ — a respectively. The 
first case is equivalent to the presence of a mass point at z(Q £ r+ , while the second 
case means that there is a mass point at — £ T~. 

In other words, any a £ D defines two orthogonal one-parameter families of 
straight lines (see figure H]): those passing through each £ £ S a in the direction 
of a/C — a, and those passing through each ( £ S a in the orthogonal direction 
— a. The points of D swept by the family {b^ ^}^es Q are the values of 6 giving 
mass points for at respectively. Hence, the values b £ D which yield mass 
points for /i ajb are those swept by {^J(es„ U {&o^}ces - Given 6, the number 
of mass points of \i a ^ is equal to the number of straight lines of both families 
{&o ^}ces a which pass through b. Moreover, each line passing through b provides 
the corresponding mass point z = =bz(£) because it crosses the unit circle at ( and 
z. This follows from 

cT (C-^) 2 IC~H 2 ,~ 

w £ T =>- — = — — -wz{Q, 



\C~a\ 



which implies that 



(C-w 



,2 



w£T\{C}, ^- -£M T w = ±z(Q. 
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Figure 4. For different choices of a, the families of straight lines 
{&ac}ces a an( j {6~^}^ gSa in blue and purple color respectively. The 
corresponding envelopes Ef (a) are the continuous curves in the same 
dark color lying on the exterior of T. The arc S a appears in grey 
color and T \ E a in yellow. The dashed curves in the interior of T 
are the envelopes Ef(a) of the families {b^ ^}c6T\s a (which are not 
depicted here). The tangent points to T split E e (a) into 3, 4 and 5 
connected components respectively from the upper to the lower figure. 
The straight lines corresponding to such components sweep different 
sectors of D. The subset of D swept by any of these sectors is the 
region of values of b giving localization for the choice of a. When 
several sectors overlap, the corresponding values of b yield as many 
mass points as overlapping sectors cover b. 
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7.2. The envelopes of {b^}^s a - The envelopes of the two families {^IceSa 
of straight lines can help us to determine the points of © swept by them. For 
the computation of the envelopes it is convenient to rewrite the equations for the 
families in a different way. Denoting A = £ — a and B = £ — b, 

^ ~ 6 ) G M T fill VWA B±V±A B±(\A\±A). 
Q-a 

Therefore, the equation for the family can be written as 

Re[(A±\A\)B] = 0. 

Remember that ( G S a is given by ( = -^e lt , cost < \a\, so A and £> can be 
considered functions of t which parametrizes the lines of the two families. Then, the 
envelope of each family is given parametrically with respect to t by the equation of 
the family together with its derivative with respect to t. This leads to the equations 

Re[X±B} = 0, Im[Y±B + X±(] = 0, (36) 

for the envelopes of 6 > respectively, where 

X ± = A±\A\, Y± = i-X ± = C±i-^-, 



only depend on a and t, but not on b. 

The system ( 1361) can be solved in B(t), thus in b(t) = ((t) — B(t), giving the 



envelopes b^(t) of the two families b^ t 



b± = ( + i Im{X± t\x ± . (37) 
a S Re(X±Y±) ± V ; 

When we let t G [0,27r], then ( runs over the whole unit circle and the two 
envelopes obviously describe a closed curve because t enters in b^(t) only 

through (. We will refer to 

E ± (a) = {b±(t):tE[0,2n]} 

as the full envelopes, to distinguish them from the original ones 

E?(a) = {b£(t):cost<\a\}, 

in which ( runs over E . The closure E e (a) = {b±(t) : cost < \a\} allows £ to 
run over the closed arc E a , that is, it only adds to Ef(a) the two limit points in 
dEf(a) = {&„(£) : cost = \a\}. Apart from being useful in some reasonings, the 
rest of the envelope E^(a) = : cost > \a\} has no interest for us because it 

comes from points ( 6 T \ S a . When referring to the set of two ± envelopes we will 
use the notation E(a) = E + (a) U E~(a), E e (a) = E+(a) U E~(a) and so forth. 

The following properties of the envelopes follow from (|37|) (see figure H]): 

• The full envelopes E ± (a) are deformed deltoides, i.e., deformed triangles 
with concave curve sides joining at three cusps. Also, E + (a) D E~(a) = 0. 

• E e (a) C C \ D, Ei(a) C 1 and dE e (a) C T. Indeed, dEf(a) = dT±. Hence, 
we will call E e (a) the exterior envelope, Ei(a) the interior envelope and 
dE e (a) the limit points of the envelope. 
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• Two contiguous cusps lie on the closed exterior envelope E e (a) if and only 
if there is a tangent point to T in the side joining such cusps. In particular, 
a cusp lies on T if and only if it is a tangent point to T, and this is also 
equivalent to stating that the cusp is a limit point of the envelope. 

Consider a given value of a. 

If the 3 cusps of one of the two full envelopes E ± (a) lie on the corresponding 
closed exterior E e (a), then E e (a) becomes tangent to T at 2 points (see the lower 
image in figure H|) . In such a case the straight lines corresponding to the points of 
Ef(a) between the two tangent points sweep the whole unit disk D and, thus, any 
value b G D gives a mass point for fi a ^ . 

In general, if T(a) are the tangent points to T of the exterior envelope E e (a), 
then E e (a) \ T(a) splits into connected components with a single cusp inside each 
component (see figure H]). The straight lines corresponding to a given connected 
component do not intersect in D, so a value b G © yields as many mass points as 
connected components have a straight line passing through b. The region of D swept 
by the straight lines corresponding to a connected component is the open sector 
limited by the tangents to the two extreme points of the connected component. If 
one of the extreme points of a connected components is tangent to T and the other 
one is a limit point of the envelope, the only straight line which limits the related 
sector of D is that one tangent to the envelope at the limit point in question. If 
both extremes of the connected component are limit points, then the related sector 
is limited by the two straight lines associated with such limit points. If a limit 
point is simultaneously a tangent point to T of the closed exterior envelope E e (a), 
then the corresponding straight line does not provide any restriction for the related 
sector of values b G D. 

In other words, each connected component of E e (a)\T(a) has an associated sector 
of values b G © to which it gives a mass point for ii a ^. The elements deciding the 
sectors are the tangent points to T of E e (a) and the straight lines corresponding to 
the limit points of the envelope, which we will call the limit lines of the envelope. 
Let us have a closer look at such elements. 

7.2.1. Limit lines of the envelopes. The limit points dE e (a) of the envelopes corre- 
spond to setting ( = = n(H ± ip a ). The related straight lines b~ ±, b + ± pass 

through such points in the directions ^C, — a and iy/C, — a respectively, which are 
the orthogonal directions \j ±ia and respectively (notice that \j ±ia points in 

the same direction as 1 ± i-rr). 

\a\ > 

More precisely, the limit lines of E+(a) are the orthogonal lines ± and the 
limit lines of E~(a) are the orthogonal lines b~ T , which are parallel to the previous 
ones. Thus the limit lines of the envelope are two pairs of parallel lines which are 
orthogonal between themselves. 

On the other hand, we know that b^ crosses the unit circle at £ and ±z(£). 
Therefore, the limit lines are: 

• +, joining to ±z((+) = ±z a respectively. 

• 6* , joining (~ to ±z((~ ) = ±z a respectively. 

7.2.2. Tangent points to T of the full envelopes. If the full envelope E ± (a) is tangent 
to T at a point £, then passes through ( in the direction orthogonal to £. 
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Figure 5. The curve in orange is the epicycloid which delimits the 
values of a with a different number of tangent points of E(a) to T. 
The cusps are at ±|. The tangents to the epicycloid passing through 
a given a cross the unit circle exactly at the tangent points of E(a) 
to T. Due to the shape of the epicycloid there are three possibilities 
according to the relative position of a and the epicycloid. These 
possibilities are shown in the figures above, where the orange points 
are the tangencies between E(a) and T. 



This means that £ _L — a or £ _L y/( — a, which is equivalent to £ || VC — a or 
C II ~ a - Thus, the tangency condition can be expressed as 

C 2 

-r 1 — € R. (38) 
C- a 

This condition is satisfied by the values of a lying on a straight line passing 
through ( in the direction ( 2 . The line picks up the values of a with E(a) having 
£ as a common tangent point to T. 

The envelope of the family of lines {a^j^eT will help us in counting the number 
of tangent points to T of E(a) for any value of a. Setting ( = e lt , t £ [0, 2ir], this 
envelope is given by equation (1381) together with its derivative with respect to t, 

Im(C-aC 2 )=0, Re(C - 2aC 2 ) = 0. 

The solution 

a{t) = (Re(aC 2 ) - Tm(aC 2 ))C 2 = (^ReC - TrnCK 2 = |e lt - \e™ 

is the envelope of {a^}^ 6 T, which is an epicycloid inscribed in the unit circle with 
two cusps at the points ±| (see figure [5]). 
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Given a £ O, the number of tangents to the epycicloid passing through a counts 
the number of points in E(a) which are tangent to T. Therefore, due to the shape 
of the epycicloid we have the following possibilities for a £ D (see figure [5]): 

(T 2 ) If a lies inside the epycicloid, E(a) has 2 tangent points to T. 

(T 3 ) If a lies on the epycicloid, E(a) has 3 tangent points to T, except at the 
cusps a — ±2, where E(a) has 2 tangent points to T. 

(T 4 ) If a lies outside the epycicloid, E(a) has 4 tangent points to T. 

7.2.3. Tangent points to T of the closed exterior envelopes. To complete the picture 
of the mass points of [i a ^ we need to know the tangent points to T of E e (a). 

We know that every line a^, £ £ T, includes all the values a £ D with £ as a 
common tangent point of E(a) to T. However, ( lies on E e (a) if and only if ( £ S a , 
i.e., Re(aC) < |a| 2 . Bearing in mind that the parametric equation of is 

a C W = C + AC 2 , A£M, 

we find that the curves separating the points of the lines {a^(A)}^ 6 T lying on the 
exterior and the interior envelopes are given by 




Re(a c (A)C) = MA)| 2 & 

The curve corresponding to A = is the unit circle, so it does not impose any 
limitation to the values a £ D. Writing ( = e lt , the remaining curve is given by 
A = — cos t, so it has the form 

a (t) = e u - e 2it cost = ^e u -^e 3l \ 

which is an epitrochoid inscribed on the unit circle with two loops and two self- 
intersections at the points ±^= (see figure [6]). 

The set A = {a(t) : t £ [0,27r]} is formed by all the values a £ D with E e (a) 
tangent to T at some point of dE e (a), i.e., at some limit point. Given a £ A, the 
number of tangencies at the limit points dE e (a) is equal to the number of lines a^, 
C £ T, passing through a, i.e., the number of tangent lines to the epitrochoid at the 
given point a. Therefore, all the values of a lying on the epitrochoid have a single 
limit point where E e (a) is tangent to T, except the self-intersections a = ±7^75, in 

which case E e {a) is tangent to T at two limit points (see figure [6] and the right 
column of figure [7} . 

For any value a £ D \ A, E e (a) has no tangent point to T on dE e (a), thus the 
tangent points to T of E e (a) must lie on E e (a). The set D \ A splits into 6 open 
connected regions. By continuity, the number of points of E e (a) which are tangent to 
T must be constant on each of these connected regions. Also, continuity arguments 
together with the fact that the envelopes E + (a) and E~(a) do not intersect, ensure 
that the distribution of tangent points between E~^(a) and E e (a) must be the same 
inside each of the above connected regions. Hence, the picture for the number of 
tangent points to T of the closed exterior envelopes E e (a) can be completed by 
simply calculating such number for one value of a in each connected region. This 
gives the following results (see figures IU E] and the left column of figure [7]): 

(Tg +1 ) If a lies inside the loops of the epitrochoid, E e (a) has 1 tangent point to T 
and it lies on one of the open envelopes Ef[a). 
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(Tg +1 ) If a lies inside the epitrochoid but outside the loops, E e (a) has 2 tangent 

points to T, one in each open envelope Ef(a). 
(Tg +2 ) If a lies outside the epitrochoid, E e (a) has 3 tangent points to T, two of 

them in one of the open envelopes Ef(a) and the other one in the remaining 

open envelope. 

When a G A we know that E e (a) has one or two limit points tangent to T, 
depending whether a ^ or a = The complete picture for a lying on the 

epitrochoid can be inferred by continuity from the previous results (see figure [6] and 
the right column of figure [7J) : 

(Tg +1 ) If a ^ ±^ lies on the loops of the epitrochoid, E e (a) has 2 tangent points 

to T, one in an open envelope Ef(a) while the other one is a limit point of 

the remaining envelope. 
(Tg +2 ) If a = ±^75, E e (a) has 3 tangent points to T, one in an open envelope Ef(a) 

while the other two are the limit points of the remaining envelope. 
(Tg +n ) If a lies on the epitrochoid but outside the loops, E e (a) has 3 tangent points 

to T, two of them in different open envelopes Ef(a) while the other one is 

a limit point. 

7.3. Localization pictures on Z + : dependence on a and b. The previous 
section shows that T 2 = (T\ +2 or T\ +2 or Tg +n ) picks up the values a G D such 

that one of the closed exterior envelopes E e (a) has 2 tangent points to T. The 
sector associated with the connected component of E e (a) \ T(a) ending in such 
tangent points fulfills D (see the cases 03, 05 and a§ in figures [6] and [7]) . Therefore, 
the region covered by T 2 , which is the closed exterior of the epitrochoid A, gives all 
the points a with n a ^ having mass points for any b G D. 

Consider now a value a G © satisfying Tg +1 . In this case E e (a) \ T(a) has 4 
connected components, each of them ending in a tangent point to T and a limit 
point. The associated sectors of D are all included in one of them (see the case a>2 in 
figures M and [7J . Therefore, such a dominant sector provides the values b G © with 
/i aj b having mass points. In the case T^ +1 we have a similar conclusion, although 
only 3 connected components appear in E e (a) \ T(a) (see the case a 4 in figures E] 
and[7J. Summarizing, the condition = (T^ +1 or Tg +1 ), which means that a lies 
on the interior of A but not on the interior of the loops, ensures that localization 
holds for all the values of b bounded by a single limit line. 

Finally, assume Tg +1 for a G D. Then, there are 3 connected components in 
E e (a) \ T(a), whose related sectors are included in that one bounded by two limit 
lines (see the case a\ in figures [6] and [7J . The points of this dominant sector are the 
values of b giving mass points for fi a ^. 

This provides the following localization criteria, which are the analogue of (A) 
and (B-t) for the case of Z + at the end of section [5] (see figures [6] and [7j): localization 
holds for the values b G D bounded by 

(Lo) No limit line <^ a lies on the closed exterior of A. 

(L x ) One limit line <^ a lies on the interior of A, but not on the interior of the 
loops. 

(L2) Two limit lines <^ a lies on the interior of the loops of A. 
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Figure 6. The left figure represents the epicycloid (in orange) and 
the epitrochoid A (in green). The self- intersections of A are at ±4g. 
For all the values of a in a given region enclosed by these cycloids the 
number of tangent points to T is constant for both E e (a) (left number) 
and Ei(a) (right number). Each crossing of the epicycloid from the 
exterior to the interior reduces the total number of tangencies to T of 
Ei(a) in 2, keeping invariant the number of tangencies of E e (a). On 
the other hand, any crossing of A from the exterior to the interior 
changes a tangency of E e (a) into a tangency of Ei(a). The number 
of tangencies to T of E e (a) is determined exclusively by A. This 
number, going from 1 to 3 as a runs over the unit disk, is indicated 
by the darkness of the green color in the right figure. 

An independent argument proves that, like in the case of Z+, for any o6D, there 
exists localization for all the values b lying on the open set S(a) defined in (1251) . 
limited by the arc S a and the straight line joining and (~ (see figure [I]). The 
reason for this is that, when £ runs over S a , the phase of £ — a performs a rotation 
of an angle it, so the orthogonal straight lines b a ^ rotate by an angle ir/2. Then, 
geometric arguments show that the pair of orthogonal families {&„^}ceSo sweep a 
region of D limited by at most two of the limit lines and that this region includes 
the set S(a). 

As a consequence, there are at most two limit lines crossing the open arc T \ E , 
and these limit lines define in any case the sector of values b G D giving localization. 
Hence, the number of limit lines crossing T \ S a is k in the case (see figures |6] 
and [7]). In any of the cases the measure fi a ^ can have 1, 2 or 3 mass points, but 
only in the cases Li and L 2 it can have no mass points. 

In the borderline case a G A the envelope E e (a) is tangent to T at a limit point, 
so that a limit line is tangent to T at such a point. Then dH a and dT a must have 
in common such a limit point (see figures [6] and [7J . 

When a is imaginary the limit lines are parallel to the real and imaginary axes 
because — a G M, and the two limit lines which are parallel to the real axis 
coincide because Im£+ = Im£~. In such the region of values of b giving 

localization becomes exactly S(a). Therefore, as in the case of Z, among the values 
of a with the same modulus the imaginary ones provide the largest region for values 
of b without localization. 
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Figure 7. Localization for one defect on Z + (from a to b). In 

red the values of b giving localization for each of the values of a in 
the right hand side of figure EJ The values of a in the first column 
are those already represented in figure HJ and the darkness of the 
color shows the number of overlapping sectors swept by the straight 
lines corresponding to the connected components of E e (a) \ T(a), so 
it indicates the number of mass points (running from 1 to 3) for the 
related value of b. The right column represents similar figures for 
different borderline situations corresponding to values of a lying on 
the epitrochoid A. The value a\ leaves two disconnected bands for b 
without localization (condition L 2 ). For a% and a 4 a unique band of 
values of b is free of localization (condition Li). The values a^, a§ and 
a 6 yield localization for any b (condition L ). Given a, the bands of 
b which are localization free are bounded by the limit lines crossing 
the arc T \ X a (in yellow). 
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8. Asymptotic return probabilities: one defect on Z + 

Like in the case of Z, we can perform the computation of the asymptotics of 
Pa s( n ) usrn g CD]) an d the canonical representative d/t = d/x a ^ instead of the measure 
dfx(z) = dfi(e~ l ®z) of the QW. From (JH]) and ( JT5l) we find that the corresponding 
functions for the state are related by tfj^^z) = t/k ^(e~ t,& z) with a = A^a 

and $ = A 2 fc+i/3. Hence, = A2fc$i*o(e~ w z), V^itX) = A 2 fc+i$iji (e~ w ;z) and 

ifTTl) can be expressed as 



^z>^( 2 )A({4)V'i,o( 



0,1^ 



(fc), 



(39) 



Again, we will omit the hat on a, f3 making the substitution a, /3 — > a, j3 at the end 
of the calculations. 



1 1 - aQ 



8.1. Masses of fx = Li a ,b- Any mass point of fx a ,b has the form 

— ± h =TT fc 1 a i A o — T— : 

Co - a 

Since h aj b(zo) = Zof aj b(zo) = 1, the corresponding mass is given by 

A(W) = v^TioT = 1 + ■ 

Performing the change of variables ((z) = —z 2 f a (z) we obtain 

C-b 
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fa.b — C " 



(40) 



(41) 



(42) 



The expression of C'( z o) remains as in (1301) for z G r+, but has opposite sign for 
z G T~ . Therefore, 



/a,&(*o) = =F 

which finally yields 

£({*o}) = - 
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(43) 
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8.2. Asymptotics of Vasi 71 ) on We need the function il> a 'p = aX + (3X 1 , 
where X (z) = 1 and X\{z) = (z~ l — b)/pb follows from the the first two equations 
oiCX(z) = zX(z). Hence, 

^*A Z ) =oc+—{z -b). 
Pb 



(0) 
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Any mass point z satisfies z f a ^ b (z ) = 1, which using PS]) gives 
Therefore, 



b=-jL. 







^i°j»fo)/Kfo})#i!o(*b) = Mi>o}) (a - Z 3 ^^! / 



^l°k^)/HW)^o°J(2o) = -£({zo}) 



« P6 I 
. a - p= — = . 



The cases with more than one mass point give in general a non convergent return 
probability p^\{ n ) due to the different factors z n appearing in fl39|) . Nevertheless, 



the case with only one mass point zq yields 

,(0) (, 



n->oc * \ |Co - b\ 



p\ 



a — (3 



Pb 



1 + 



Pi 



1 + 2- 
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Co-b 



Co-b 

2 



Then, all the states at the origin exhibit localization except that one defined by 



j3 = a- 
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